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ABSTRACT 

We show how the excursion set moving barrier model for halo abundances may be 
generalized to the local non-Gaussian /„; model. Our estimate assumes that the dis- 
tribution of step sizes depends on but that they are otherwise uncorrelated. Our 
analysis is consistent with previous results for the case of a constant barrier, and high- 
lights some implicit assumptions. It also clarifies the basis of an approximate analytic 
solution to the moving barrier problem in the Gaussian case, and shows how it might 
be improved. 
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1 INTRODUCTION 

Detections of non-gaussianity can discriminate between different inflation models (e.g. Maldacenal 20031 ) . The local /„; model, 
where the primordial perturbation potential is 

$ = <^+/„,(<^'-{<^')), (1) 
where is a Gaussian potential field and fni is a sc alar, has been the subject of much recent study fe.g. jBuchbinder et al 



Khourv fc Piazza' 2008; Silvestri & Trodden '2OO8, and references therein). Constraints o n this mo del tend to be of two types 



from the CMB (Komatsu _et^_aL 



20O9j_^Hikagc ct_al. 



2008: Yadav fc Wandch 2008.: McEwen et al.l,2008.) and from large scale 



structures in the Universe jKovama et al] 19991: Matarrese et al. 2000l: Iscoccimarro et al.ll20oJ 



Sefusatti fc Komatsu 



200i 



Izumi fc Soda"2007'; 'Lo Verde et al."2008'; 'Dalai et al."2008'; 'Matarrese fc Verde"2008'; 'Carbone ot al, 
2008 : Slosar et al. .200 8; McDonald 2008; Taruva et al. 2008 : Slosa r 2008; Grossi et al. 200 8: Kamionkowski et al 



2007 



'8': ' Afsh ordi fc ToUev 



Desiacaues et ahlboosl : iLam fc Shethli2009l : Icrossi et al.ll2009l : I Valagea3 120091 ) 



200S 



One of the fundamental quantities of interest in such studies is the abundance of virialized dark matter halos. lPress fc Schechter 



1 19741 ) suggested that the abundance of collapsed virialized halos may be estimated from the statistics of the initial fluctu- 
ation field. They used the assumption that halos form from a spherical collapse to argue that such o b jects started out as 
sufficiently overdense regions in the initial fluctuation field. The excursion set approach of iBond et al. ( 1991 ) allows one to 
estimate halo abundances in Gaussian theories; in this context, the spherical collapse model is associated a barrier of constant 



height. If the collapse is triaxial, then the barrier height is stochastic with a mean that is not constant ijSheth et al.l 12001 



Ignoring the stochasticity but including the changing of the barrier height allows the excursion set approach to provide a 
simple parametrization of the effects of triaxial collapse on halo abunda nces (Sheth fc Tor men 2002! ) . The main goal of the 
present work is to show how to generalize the moving barrier formulae of Isheth fc Torrnen (2002) to the local non-Gaussian 
/„; model. 

Section [2] provides explicit expressions for the one- and two-point distribution of the overdensity in /„; models. These 
are used, in Section [3l to estimate how the mass function of virialized objects is modified when /„; 7^ 0. This section also 
clarifies earlier work on the Gaussian (/„; — 0) case. A final section summarizes our results. Two Appendices provide useful 
approximations and other technical details. 



* E-mail:tylam@sas. upenn.edu, shethrk@physics.upenn.edu 
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2 THE LOCAL NON-GAUSSIAN MODEL 

We are interested in models where the primordial perturbation potential is given by equation ([T}. We will use -P^(fc) to 
represent the power spectrum of (f>; in what follows we will set P4>{k) = Afc"^"'', where Us ~ 1, and A is a. normalization 
constant that is fixed by requiring that the rms fluctuation in the associated non-Gaussian initial density field (which we will 
define shortly) has value erg. The power spectrum and bispectrum of the "1> field are 

P*(fc) = P^{k) + ^ J'dq [PMPA\k -q\)- PAk)PM " P^{k)Pd\k - q\)] , (2) 
B#(fci,fc2,fci2) = 2Ui [P^(fci)P^(fc2) + cyclic] + 0(7^0 (3) 



i Scoccimarro et al.ll2004l ) 



Most of the complication in /„; models arises from the fact that we are almost always interested in spatially smoothed 
quantities. In particular, the quantity aSa = {S'^\R) / {5'^\R)^^^ , will play an important role, because it represents the leading 
order contribution to the non-Gaussianity (note that it is proportional to fni)- Fortunately, smoothing is a linear operation, so 
the smoothed variables are just linear combinations of the unsmoothed ones. Hence, if W{kR) denotes the Fourier transform 
of the smoothing window of scale R, then, to second order in /„;, 

{5^\R) = a\R) ^ f ^ATTk''M^{k)P^{k)W\kR), and (4) 
(27r)'^ J k 

{5'\R) = 2/„, / ^ktM{ki)W{kiR) f ^klM{k2)W{k2R) / d^ia fcfa M(fci2)VK(fci2P) -^"^^1' ^'"^ (5) 
(27r)4 J ki J k2 J 2f„i 

where M{k) = {'iD{z)(?) / {^Q,mH'o)T{k), T{k) is the CDM transfer function and D{z) is the linear growth function. In 
hierarchical models, a and aSa are both monotonically decreasing functions of 7?; this will be important in what follows. 
Appendix |X] provides a useful fitting formula for a Sa, and shows that it is only a weak function of scale. 



2.1 Edgeworth approximations for p{5\R) and p{5, A\r, R) 



Because we are interested in small departures from Gaussianity, the Edgeworth expansion provides a convenient form for the 
distribution of S smoothed on scale R: 



p{5\R) dS : 



a{R) 



2tt(t{R) 



■dS - 



Po(5\R) dS, 



(6) 



where a{R) is given by equation Q, aSa = {x^)/{x^)^^^ = 2/„;7^/cr^, and H3{u) = u{u^ — 3) with u = 5/a{R). Because H3 
changes sign sd, u = \/3, the Edgeworth expansion is not always positive, making it ill-suited for studying (at least one of) the 
tails of the distribution. The expansion becomes negative when {aS3/Q)u{h'^ — 3) = —1, meaning !/(i^^ — 3) = —200 (O.OS/aSs), 
and note that aSs has the opposite sign to /„;. For /„; ~ 100, we have v{u^ — 3) = 200, so there certainly are problems at 
v > 6, making the expansion suspect at slightly smaller values. For /„; ~ —100, there are problems at i/ < —6. Because the 
Gaussian piece falls exponentially with , the fact that the expansion may not be accurate at large u may not matter - but 
for l arger \fni\ values, th is limitation of the Edgeworth approach should be borne in mind. 



Lam fc Sheth 1I2OO9I ) have used this expansion in the context of modelling the one-point distributi on of the evolve d 



nonlinear 5, where a S3 is smaller, and large u values are indeed rare. Here however, we will follow Lo Verde et al. ( 20081 ) . 
and use it in our model of halo abundances. In this case, it is the large 1/ tail which is of most interest, and this is precisely 
where the Edgeworth expansion is most suspect. However, note that, when modeling halos, one is most interested in the 
regime where > 0, so, for /„; < 0, the Edgeworth expansion is positive definite except for large a S3, for which we know 
the Edgeworth expansion is not useful anyway. In any case. Appendix |A] shows that large values of aSs are not currently a 
concern. 

For reasons that will become clear shortly, we will also be interested in the value of the field when it is smoothed on two 
different scales. For small departures from Gaussianity, the bivariate Edgeworth expansion should provide a good description. 
It is 



p(p,u) = 

i Kotz et al.ll2000l ). where 

^ 5 

^= (52)1/2' 



1 



27ryr 



: exp 



2(1 -g2) 



1 



+ 



^2lH21 + Ai2i/l2 



(7) 



(A2)l/2 



q = {fiu), A„„ = (^""i^") c, and Hmn{^;y,q) = 



(1-^2)2 ' 
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with 

^30 (/i, ^,<l) 



(u — qu)^ 

ho3{v,^.,q) = — 5 3(^i - qv), 

i — 



h2i{fJ,,iy,q) = hu{v,^J.,q) = 2q{^- qu) - [v - q^) + 



{v-q^){^-qvf 

1-02 



If our convention is that A is the field on the larger smoothing scale, then, to lowest order in A03 = ctSs, 



exp 



{^i-qvf 



2(1 



1 + 



6 



6 



(8) 



In what follows, we will set q^ — {A'^)/(5^) (as it is for a Gaussian field smoothed with a tophat in fc-space; this is standard 
for the excursion set approach). 



3 HALO ABUNDANCES 

Recall that ct is a monotonically decreasing function of smoothing scale R. In the initial conditions, where fluctuations are 
negligible, the smoothing scale R contains mass m = pAivR^/S almost surely, so a, R and m are equivalent variables. Now, 
let dn/dm denote the comoving number density of halos of mass m, and let p denote the comoving density of the background. 
If /(m) denotes the mass fraction in halos of mass m, then 

F(<a(R^i3m/4npr^'))=Fi>m)= H dm fim) = H^^, (9) 

An estimate of F, then, is an estimate of dn/dlnm. The following sections describe the excursion set estimate of F for 
Gaussian initial conditions, and how this estimate can be extended to the /„; 7^ models. 



3.1 Excursion set approach 

The excursion set approach ( Bond et al. 1991 : Lacev fc Cole 19931 : Sheth 1998l l relates the number of haloes of mass m to 
the first crossing of a suitably chosen barrier b{a) by a suitably chosen set of walks. The simplest implementations of this 
approach consider the first crossing of 6((j) by an ensemble of uncorrelated random walks with uncorrelated steps. While 
neglecting both types of correlations is far from ideal, previous work shows that this allows one to write down simple analytic 
expressions for the first crossing distribution and how this distribution is related to halo abundances, which together provide 
reasonably accurate descriptions of halo abundances as well as their formation histories. 

The barrier shape 6(a) is set by the physics of collapse. The spherical collapse model has a barrier of constant height 
b((T) = 5sc, whereas barriers of the form 



(10) 



with — 0.4 and 7 — 0.6 may be related to models in which halos from from a triaxial collapse ijSheth et al.l 12001 



Sheth fc TormerJ 120021 ). The physics of collapse has a — 1, but setting a = 0.7 results in a predicted dn/d\nm which is 



in much better agreement with the abundance observed in simulations. (See lSheth et al.ll200ll . for further discussion of why 
a /I). 



3.2 The Gaussian case: Spherical collapse 

Let p{5, s) denote the probability that a randomly placed cell in the initial distribution has overdensity S when the smoothing 
scale is such that {5^) = = s. Classify all cells by the largest scale 5* < s on which they had overdensity greater than Sc. 
Then, provided 5 > 5c, 

p{S, s) = f° dS f{S, 5c) p{5, s\5c, S, first), (11) 
Jo 

where / is the fraction of cells for which S was the largest smoothing scale on which the overdensity was greater than Sc, 
and p{S,s\5c, S, Hist) is the probability that the overdensity on scale s is 5 given that 5* was the largest scale on which the 
overdensity exceeded Sc. If one views a plot of S versus smoothing scale as something which resembles a random walk, with 
large smoothing scales to the left, then S is the first 'time' that the walk crosses Sc'. hence the word 'first' in the expression 
above. Therefore 

poo pS pOG pS 

P{Sc,s)= dSp{S,s)^ dSf{S,Sc) dSp{S,s\Sc,S,Rrst) = dS fiS,Sc) PiSc,s\Sc,S,&rst). (12) 
Jsc Jo Jsc Jo 
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In what follows, we will use the subscript '0' to denote quantities associated with Gaussian initial conditions (for which 
fni ~ 0). For a Gaussian field smoothed with a tophat in fc-space, 

Po('5i,sl(52,S',first) ^ po{5i,s\52,S) = poiSi -62,8- S), (13) 

so Po{5c, s\5c, Sj&Tst) — Po(0, s — S) = 1/2, and equation ([12} implies 

Po(.., .) . f dSMS, s)^iysM^.^^ (14) 
Differentiating both sides with respect to s shows that /o is simply related to po: 

dPo _ f d{5^/VJ}_ £^_Ms^ 



3.3 The Gaussian case: Ellipsoidal collapse 



If Sc depends on s, as in some parametr i zation s of triaxial collapse, then this simplicity is lost. In particular, for the barrier 
given in equation (fTO)). Sheth fc Tormen 1 2002^ show that 



sfo(s,b) ^ ^ — exp 



/27rs 



"2s 



1 + 0.067 



2^7 



To see where this comes from, note that the analogue of equation (|14p becomes 



ds 



-b''/2s 

PhTs 



db b_ 

as " 2s 



Ms,b) 



dSMS,B) 



-(b-B f /2(s-S) 

^/Ms - s) 



db [b - B) 



ds 2{s - S) 



(16) 



(17) 



However, the first term in square brackets on the right hand side is independent of S, so this term integrates to p{b, s) db/ds. 
Since there is a similar term on the left hand side, this leaves 



b e-"'/^- 

S J2-KS 



= Ms,b)+ / dSfo(S,B) 



-{b-Bf/2{a-S) 



V27t{s ~ S) (s~S) 



where we have multiplied both sides by a factor of 2. Writing -8(5*) as a Taylor series around b{s) implies 



be-"'^'^ ,,^^d^b 
= /o(s,6) + }^ — 



/27rs 



dSfo{S,B). 



-{b-BY /2{s-S) 



\/27r(s - S) 



(18) 



(19) 



sfo{s,b) 



db 

OS 



/2^ 



-{b-By /2(s-S) 



VMs - S) 



{S/s - 1)^ 



(20) 



Equation (|16|l. the approximation of Sheth fc Tormei] l|2002l ). corresponds to ignoring all S/s terms and then keeping only 
the first few terms in the series. 



3.4 Dependence on /„( 

We now extend the analysis above to the case in which /„; 7^ 0. Our primary assumption is that the expressions above remain 
valid if p on the left hand side of equation p2|) is given by the Edgeworth expansion, and the conditional distribution on 
the right hand side is given by the bivariate Edgeworth expansion. This is not quite right, since the steps in the walk are 
now correlated, and we are assuming that the statistics of walking from {B,S) to (&, s) do not depend on the fact that the 
walk did not cross B before 5*. In efi'ect, this means that we assume that the extra constraint that S' < 5c for all S' < S 
makes no difference, so the only difference from the Gaussian calculation is that all PDFs are now replaced by the appropriate 
Edgeworth approximation. 

To see what this implies, it is convenient to first define 



in = d5po{S + b,s\B,S)hm.n ^ 
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where 



Go3 



1 - 



Po 



B 



H3{B/VS)Po 



b-B 



1 



1 - 



{S{s - S){2S - 3s) + [5^(6^ + bB + B 

b-B 



3sS{bB + B^)+3B^s'^]}po 



B 



B? B{b-B) 

-S 



Ql2 ~ 



-sS{s-S) + {bS-Bs)' 



■po 



S 
b-B 



Po 



In addition, define 



POO 

Ga = / dSpoiS + b, s\B, S)//3(B/%/S) = Ha{B/Vs)Po 
Jo 



b-B 



If we ignore the scale dependence of a S3 (Figure [T] shows that it is weak), then the derivative of equation p2p w.r.t. s is 



(22) 
(23) 



dPojb/^) 

ds 
where 



H3 



.ns,b) 



-f 

Jo 



dS/(S,B)^Po 



b-B 



+ 



6 Jo 



dsf{s,B: 



d_ 

ds 



£{s,S)po 



b-B 



£{s,S) = 



B^ + bVs) 



(b- g)'v/s~g 



(24) 



(25) 



Notice that, in contrast to the case in which /„; = 0, the final term on the right hand side of equation (|24|) is non-trivial 
because it is now s-dependent. 

Substituting the /„; — solution for the pieces with subscript zero yields 



Ms,b) 



1 + 



6 



f{s,b) 



1 + 2 



+2 



d 

dS—Po 

OS 



b-B \ fiS, B) - MS, B)[l + {aS3/6)H3{b/^s)] 



f{s,b) ds 



fis,b) 



£{s, S)po 



B 



(26) 



This is an integral equation for f{s, b) which is valid when fni 7^ 0. Clearly, the zeroth-order solution is simply the f„i — 
solution (the Gaussian case) times a correction term which depends on a S3 and on the barrier shape b{s). We can include the 
next-to-leading order contribution as 



/(s,&) = /(«)(s,6)(l + 



f(0)J 



Ms,b) 



6 



_b_ 



6 



G{s,b) 



(27) 



where fo is the first crossing probability associated with uncorrelated steps when /„; — (approximated by equation ll6p . 



gis,b) = 2 dS 



fo[s,b) [ds [ \^s- S 

where only terms to first order in crSa have been kept, and 



OS 



B 



3B 



d£_ 

9s ^^i~5(V^+^)2 ' SV^~^{^s + ^) 



9B 



B^s + by/S 



3B{y^ - VS) 

2^s 

b-B 



+ {B + b) 



5 



[H3{B/Vs)-H3{b/V^) 



+ Vs^iV~s-^) 

OS 



Thus, 



f{s,b)^fois,b) 



'2(s-S)^ 

OS 



0-S3 



b-B 



3 (b-ByvT^ 

2 ^{^+vsr ■ 



(28) 



(29) 



(30) 



The analysis simplifies somewhat for a constant barrier - see Appendix [B] Previous a nalyses of the case in w hich the barrier 
is a constant have either explicit ly ignored the final te rm in the square brackets above ( Matarrese et al. 2000l ) ■ or have missed 
the fact that it is there entirely i Lo Verde et al. 20081 ). 

Figure[T]shows the effect of Q{s, b) on the ratio of the first crossing distribution when /„; — —100 to that when /„; — 0. This 
is done for three different barriers: the top curves (offset upwards by 0.4) show results for the moving barrier of equation (|10p . 
the middle are for a constant barrier of height b = Sc, and the lowe st (offset downwards by 0.4) for b = y/0.7Sc . The middle 
curves are for the barrier shape that has been previously studied I Matarrese et al. 2000l : Lo Verde et al. 20081 ) : differences 
between these and the bottom curves are due to the barrier height. Differences between the top and bottom sets of curves 
are due to the s dependence of the barrier. 
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Figure 1. Ratio of the first crossing probability for /„; = —100 to that when /„; = 0. SoUd (black) curves set aSs = 0.022 and 
Q{s,b) = 0; dotted (red) curves include Q{s,b) but Iteep aSs fixed, and the dashed (green) curves include the scale dependence of 
aSs but set Q(s,b) = 0. The top set of curves (offset upwards by 0.4) show results for a moving barrier (equation HOC : middle curves use 
a constant barrier with Sc = 1.66 and lowest curves (offset downwards by 0.4) use a constant barrier of height \/0.7<5c. 



In each case, the black (sohd) curves show the resuh of neglecting the factor of G{s) and setting a S3 = 0.022 (i.e. we ignore 
the fact that it depends weakly on scale). The red (dotted) curves include the fact that G{s, &) 7^ (or the corresponding G(s) 
(equation IB8[) for the constant barrier) with the same (fixed) value of a S3. The differences between these two sets of curves 
are small for all three barriers, indicating that the contribution from the factor Q is small. The effect of again ignoring Q, but 
now including the scale dependence of a S3 (using equation lAlfl , but setting da S3 /ds « is shown by the green (dashed) 
curves. 

Figure [2] shows the ratios of the first crossing distribution for different barrier shapes. The scale dependence of aSs is 
included and computed by the approximation formula (equation IA1[) . Black (solid), green (dashed), and cyan (dotted) curves 
plot the ratios for moving barrier b{s), b — Sc, and b — V 0.7 Sc respectively (all without Q). The red (dot-dashed) curve shows 
the term in square brackets in equation (I30|l : this includes the scale dependence in the factor a S3 but it sets dlnaSs/ds ~ 0, 
and it ignores the fact that this scale dependence will also modify the Q{s,b) term. The inclusion of the factor Q(s,b) has 
a small effect compared to using different barrier shapes. The three barrier shapes result in slightly different predictions for 
how fni modifies halo abundances. Therefore, comparisons with measurements in numerical simulations may indicate which 
barrier shapes better describe halo formation. 

Figures [1] and [2] show that the effect of G{s, b) is small compared to the effects of including the scale dependence of aSs, 
and the effect of using different barrier shapes. Therefore, it is a good approximation to set 

b 



fis,b)^Ms,b) 



1+^H3 
D 



(31) 



When the barrier is constant, b = Sc, then this reduces to the expression presented by Lo Verde et al. 1 2008h . Our analysis 
shows that there is an additional correction factor which their derivation missed (our factor of (crS'a/G) G), but that this 
happens to be small. On the other hand, they include a term which comes from the scale dependence of d{aS3/6)/dlns which 
we are ignoring (we have checked that it is small). 
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Figure 2. Ratios of the first crossing probability for /„; = —100 to that when /„; = 0. Black (solid), green (dashed), and cyan 
(dotted) curves plot the ratios for moving barrier b{s), b = Sc and b = \/ 0.7 5c when G{s) is set to 0. The red (dot-dashed) curve shows 
equation II30I I. 



4 DISCUSSION 

We showed how the excursion set approach may be extended to model halo abundances when the initial conditions were 
non-Gaussian. In this approach, the estimate of halo abundances is related to the first crossing distribution of a suitably 
chosen barrier by a suitably chosen set of walks. The physics of collapse sets the barrier shape, and the statistics of the 
initial fluctuation field sets the properties of the ensemble of walks: e.g., how steps in a given walk are correlated, and 
whether the walks are independent of one-another. Our analysis assumes that the steps in a walk are uncorrelated, and the 
appropriate ensemble contains all possible independent walks. The first assumption may be well-motivated only for Gaussian 
random fields - we argue that it may remain a useful approximation for weakly non-Gaussian fields. The second assumption 
ignores the fact that averaging ove r the full ens e mble of uncorrelated walks is only an approximation to the more physically 



appropriate ensemble described in ISheth et al.l (|2001h . Previous work suggests that this allows one to write down simple 
analytic expressions which are reasonably accurate. 

Our analysis was done in two steps - the first showed how the calculation depends on the physics of collapse: spherical 
and triaxial collapse models are associated with 'constant' and 'moving' barriers (Section [3]). For Gaussian initial conditions, 
our results clarify the nature of approximations made in previous studies of the moving barrier model (see Section I3.3|l . 

For non-Gaussian initial conditions, our analysis assumes that the correlated nature of the steps (in non-Gaussian models) 
changes the step-size distribution in a calculable way, but that steps are otherwise independent. For weak non-Gaussianity, 
we approximate the change to the step-sizes by using the Edgeworth expansion (equations |6] and (Tjl This is not strictly correct 
- it is an approximation which may be accurate for weakly non-Gaussian fields. (In this context, our analysis of the constant 
barrier model showed that previous work on this problem, which made the sa me assumptions as we do, had missed some 
terms. However, these turn out to be small.) As we were completing this work, iMaggiore fc Riottd (2009|) presented a very 



different analysis of the constant barrier problem which yields consistent results, suggesting that our neglect of the additional 
correlations between steps associated with /„; 7^ is reasonable. 
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Ours is the first analysis of tlie moving barrier problem for non-Gaussian models: the moving barrier yields quantitatively 
different predictions for halo abundances th an does the const ant barrier model, at a level that current simulations should 
be able to detect. This is explored further in Lam et al. ( 20091 ) . So we hope that our results will benefit problems which use 
halo abundances to constrain the natur e of the initial fluctuat ion field. They also provide a key ingredient to the Halo Model 
interpretations of how galaxies cluster (jCoorav fc Shethlbood ). 
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Figure 1. The quantity aS-j,{r^ for /„; = 100. Dott ed (black) curve is the nu merical result, dot-dashed (cyan) curve is equation I IAlll . 
and dashed (green) curve is the approximation from lScoccimarro et al ] l|2004l ). which is only expected to be valid on large scales. 



APPENDIX A: APPROXIMATION FORMULA FOR oSz 

The quantity aS-i measures the strength of non-Gaussianity in the smoothed field. Unfortunately, it requires the computation 
of several numerical integrals. We have found that the following provides a good approximation: 



o-Ss ~ 12 /„ 
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(Al) 



where the first term on the right hand side is the large scale approximation given in Scoccimarro et aP 1 2004h . If we set 
ro = O.lMpc/ft, mi = 0.5, m2 = 0.7, then (jS3.{ro) = —0.0316, and equation (|Aip is rather accurate (see Figure[l|. 



APPENDIX B: HALO ABUNDANCES IN THE CONSTANT BARRIER MODEL WHEN Fnl / 

The main text considered the general case of a moving barrier. When the barrier is constant then it is convenient to define 
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If the scale dependence of cr53 can be ignored, then A„ 



■ <t53, so, to first order in a53. 



p{5 + 5c,s\S,,S) = ^ + -^^Eis,S), 
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Note that E{s, 5) = when s = 5, so when cr53 is a constant, then 
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The main text assumes that the second term on the right-hand side is neghgible compared to the first. In the approximation 

where 0-53 is constant, this makes f{s,5c) = /()(s,(5c) [1 + {crSs/Q) H3{Sc/(j)]. To sec if this is accurate, one can substitute 
this expression for f{s, 5c) into the integral, and check that this contribution really is negligible compared to the first term, 
/(5,5c)/2. Namely, write 
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then note that, to leading order in ct53, wc may approximate /(5, 5c) ~ /o(5, 5c), so the second term in square brackets is 
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